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I. INTRODUCTION 



A. Background 



Recent studies 12, y, |4( have reported observations of axial segregation of particles in a 
suspension contained in a horizontal cylinder rotating about its axis. While the segregation 
of particles in flowing suspensions has been widelystudied in a variety of flows, there is an 
important difference in the geometry studied by Q, 3, 0, Q], namely, the presence of a free 
interface of the suspension: segregation was found to occur only when the cylinder is par- 
tially filled with suspension. These experiments have led to attempts at theoretical analysis 
Q], further experiments Jf|, and criticism [7]. As in earlier experiments on shear-induced 
migration, the suspension comprised small neutrally buoyant particles dispersed in fluid of 
viscosity much higher than that of water, so that the Reynolds number (based on particle 
size) was negligibly small, and the Peclet number was very large. The former implies the 
irrelevance of particle and fluid inertia in determining the state of the system, and the latter 
means that Brownian motion of the particles is negligible. Particle segregation in suspen- 
sions has important consequences in industrial processing, where it is usually important to 
disperse the particles uniformly to obtain a product of predictable quality. It is also possible 
that segregation could be used to separate particles from the fluid, or particles of different 
sizes from each other. From a fundamental perspective, it is important to understand the 
physical or mechanical origin of segregation. 

In a recent paper j^j], we had proposed a model to explain the observations in [li 13, Q ; 1^] based 
on a simple and intuitive picture of the dynamics of the suspension film coating the wall of 
the cylinder. We proposed that a balance between gravity, viscous shear forces, and surface 
tension determine the evolution of the film thickness, which in turn determines the shear rate 
in the film. The formulation was completed by coupling this balance with the conservation 
equation for particle mass, in which the particle concentration evolves according to shear- 
induced diffusion jsj]. Our analysis was simplified by the assumption that the variations of 
the flow variables within the cross section (i.e. in the r-9 plane, see Fig. H} at any given 



2 



axial coordinate are unimportant, and that segregation is driven by axial perturbations of 
the particle concentration and the film thickness. We conducted a linear stability analysis, 
and predicted that the homogeneous suspension is unstable to small perturbations provided 
a certain discriminant, determined by the dimensionless shear- induced diffusivities and the 
effective viscosity of the suspension, is negative. An important feature of our model is the 
necessity of a free interface, which is in agreement with the observation of 
One of the aims of the present paper is to lay to rest once and for all the doubts raised 
- in an ad hominem comment \t\\ - about the validity of our approach. From a complete 
2D analysis, we are able to show where Eq. (4) of p, which is the force balance mentioned 
above, comes from. Incidentally, this part of the derivation was already presented in our 
rebuttal to the comment. 



B. Summary of results 

Before describing the model and our calculations in detail, let us summarize the achievements 
of this paper. We consider the coating of a plate, with coordinates in the (x,z) plane, being 
pulled out of a pool of suspension at constant velocity in the x direction (see Fig. EJ). From 
he Stokes equations and the particle mass balance, as determined by shear-induced diffusion 
11, we derive a pair of coupled partial differential equations for the film thickness and the 
particle concentration as functions of time t and the spatial coordinates x, z. From these we 
obtain the steady-state thickness and concentration profiles (see Fig. EJ) and the equations 
for the evolution of small perturbations about this steady state (see Eqs IHB^l . The form 
of these equations, for variation solely along z, and ignoring surface tension, is identical 



0. 



Two instability 



to that presented on general physical grounds in |5|, and criticized in 
mechanisms emerge from this analysis - the axial instability proposed in our earlier one- 
dimensional treatment jjj], and another "pile- up" growth along the x direction (see §IIIj) . 
As stated in our earlier work [5,l], we find that if the discriminant E, which depends on 
the dimensionless shear induced diffusivities, is below a critical value E c , the growth rate 
varies as k 2 z for small wavenumber k z in the z direction, turning over at larger wavenumbers 
on a scale determined by surface tension. There are two differences between \^ and the 
present results: (i) Because the eigenfunctions of the growing modes all have significant 
x-dependence, the circumferential instability always acts, so that the growth rate is nonzero 



3 



FIG. 1: Schematic diagram of film coating in a horizontal rotating cylinder 



even for k z = 0. The axial instability of rides on this background "pile-up" growth-rate. 
As a result, we predict E c > whereas the one-dimensional analysis of \^ predicted E c — 0. 
(ii) The wavenumber g* of peak growth rate varies as the inverse fourth root, not square root, 
of the surface tension. Figs 0] and 03 show the effects of changing the surface tension and the 
discriminant, respectively, on the instability growth rate as a function of k z . These minor 
differences apart, the mechanism proposed in |5[ for the segregation instability of Quag 
remains viable. 

The remainder of this paper is organized as follows. Starting from the Stokes equations 
and the kinematic condition at a free interface, we first derive in [(n] the equation for the 
dynamics of the film thickness, making only the lubrication approximation as in almost all 
studies of the drainage of viscous films. We show that this equation reduces to Eq. (4) of 
(save the surface tension term) if certain simplifying assumptions are made on the variation 
of the film thickness in the direction of wall motion. In §1111 we conduct a linear stability 
analysis of this equation, coupled with the conservation equation for particle mass, and show 
that it yields essentially the same result as that the simple model we had proposed in P|. 
We close the paper with a summary and conclusions in £UVI 



II. DYNAMICS OF A SUSPENSION FILM WITH A FREE INTERFACE 



We consider the coating of a flat plate rising 
the vertical, shown in Fig. EJ Several studies 



±rom a poo 

3, In, 



12 



of suspension at an angle 9 from 
131 ] have investigated the classical 



problem of coating of a flat vertical plate rising from a pool of liquid. The object of all 
these studies was to determine the asymptotic film thickness far from the pool, as that is 
the primary variable of interest in industrial coating operations. The dynamics and stability 
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FIG. 2: Coating of a vertical plate rising with velocity v w from a pool of suspension. 

of coating flows of pure liquids also have been analyzed 1171 1 1811 . in an effort to explain 
instabilities and pattern formation that are observed |igL I20L I21L 22 1. In this study, we 
consider the dynamics of the coating flow of a suspension whose viscosity is a function 
of the particle concentration. Therefore, the dynamics of the film thickness is coupled 
to the dynamics of the particle concentration. The particles are assumed to be neutrally 
buoyant, i.e. of the same density as the liquid. The pool is assumed to be large enough that 
the particle volume fraction concentration remains constant at <p p . We ignore variations 
in particle concentration across the thickness of the suspension film, and assume <fi to be 
independent of y. 

Though our intention is to address the observations of Q, 0, 0, 0] in rotating horizontal 
cylinders partially filled with suspension, we choose to study the problem of a flat plate 
rising from a pool, as it simplifies the analysis. While the equations describing the dynamics 
of the film in a curved wall which joins the suspension pool at the far end are almost identical 
to the ones we derive here, their numerical solution is a greater challenge. It is clear, however, 
that the two problems are qualitatively similar, and are identical in the limit of the cylinder 
radius becoming very large compared to the film thickness, provided the fraction of cylinder 
filled with suspension is not small. 

From earlier investigations on coating flows, 10, 11, 12, it is clear that the liquid film 



coating the plate can be divided into three regions (see Fig. |2J): (a) the static region at the 
bottom, where there is a balance between surface tension and gravity, (b) the "entrainment" 
region where surface tension, gravity and the viscous shear stress all play a role and (c) the 
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region of constant film thickness at the top where gravity is balanced by the shear stress at 
the wall. The solution of the problem is obtained by matching the solution of the entrainment 
region wit 
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11 



12 



i that of the static region at one end and the straight film region at the other 
This solution yields, among other things, the asymptotic film thickness w c . 
It is in general a function of the Capillary number Ca = v w r//a; it varies as Ca 1//6 when 
Ca <S 1, and becomes independent of it when Ca 3> 1 



In the static region, the balance of surface tension and gravity yields 



w 



.'2 



w ■ 



3/2 



pg cos 9, 



(1) 



the primes indicating differentiation with respect to x. Integrating the above with the 
boundary condition u>' — > — oo at x = j| gives 

w' pg cos 9x 2 



w 



12 



1/2 



2(7 



(2) 



The point at which the slope vanishes is x = {2a / pg cos 9) 1 1 2 . Following the prescription of 
jlfl| . we match the curvature at this point, 
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2pg cos 9 



1/2 



a 



(3) 



with the curvature at the lower boundary of the entrainment region. 

In the region of constant film thickness, w = w c , but its value is not known a priori - every 
value of w c satisfies the balance between gravity and the wall shear stress, each corresponding 
to a particular value of the mass flux. The value of w c must therefore be obtained by 
matching with the solution in the entrainment region. Thus, the static and constant film 
thickness regions serve to set the boundary conditions for the entrainment region, and it 
is the solution in the entrainment region that determines the asymptotic film thickness w c 
U> 3 [3 and, as we shall see in what follows, the film dynamics. 



We now consider the entrainment region. As stated earlier, we are interested not just in the 



ows variations in the 
12 , [k| , we simplify 



steady state, but in obtaining a dynamical equation for w that also al 
z direction. Following earlier studies on coating of a flat plate [1 1 
the momentum balances by invoking the lubrication approximation, on the assumption that 
the film thickness w is small compared to the distance in the x direction over which there is 



variation of w. The x- momentum balance for the suspension in the film reads 

d 2 v x dp 
r] s — — = — + pgcos9. 
ay ox 



(4) 



Integrating the above with the no-slip boundary condition at the cylinder wall (y = 0) and 
the stress free boundary condition at the liquid-air interface (y = w) gives 



1 dp y 2 

— (pg cos 9 + ^-)(— - yw) + v w , 



Vs 



dx jy 2 



(5) 



where v w is the velocity of the plate. The continuity of normal stress across the interface 
yields the relation 



dv„ 



,d 2 w d 2 w. 
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(6) 



dy v dx 2 ' dz 2 

where we have again invoked the assumption that the gradient of w is small, i.e. 
\dw/dx\, \dw/dz\ <C 1, in linearizing the expression for the curvature of the interface. The 
contribution of the viscous normal stress r] s dv y /dy was ignored by l| Q], but retained by 
Q. Retaining it improves the agreement of the predicted asymptotie eoating th.chness with 
experimental observations We drop this contribution, thus considerably simplifying the 
analysis; if it is retained, the pressure cannot be eliminated from the equations of motion, 
and hence an additional field must be solved for. Substituting © in (JTj) with the above 
simplification yields 

pg cos 9 a ( d 3 w ( d 3 w 

v 



' 2 



yw) + v w . 



(7) 



rj s T] s \ dx 3 dz 2 dx 
In the z- direction, flow is driven only by surface tension due to a gradient in the curvature. 
The balance of momentum for this direction reads 



d 2 v 7 _ 



u v z i d w crw^\ , . 

Vs W = ~ a \dx T d~z + ~d?)' () 

which upon integration, again using the no-slip boundary condition at the wall and no shear 
stress condition at the interface, yields 
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v z = -- 



d 3 w d 3 w s 
+ 
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- yw). 



rj s \dx dz dz 6 J y 2 
The velocity v y is determined from the incompressibility condition, 



dv„ dv 7 



(9) 



(10) 



dy dx dz 

Upon substituting ((7|) and in (fTU|) and integrating, using the no-penetration boundary 
condition at the wall, we get 

y 2 dw a 
~2~dx 



pg cos 9 a id w d w 
dx 3 dz 2 dx 
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dx 2 dz dz 3 



1 6 



y 2 w. 



(ii) 

(12) 



Finally, the kinematic constraint on the interface that its velocity be equal to that of the 
suspension adjacent to it yields 



dw dw dw 

■dt +V 'te +V 'te =V * aty = W - (13) 



Substituting (JZj), © and JEJ) in (JTHJ) , we have 



dw dw d 
dt w dx dx 



d 3 w dV 



pg cos 9 — a I 2 — h 



dz dx dx 



vf \ _ d_ 
3r] s J dz 



aw 3 ( d 3 w d 3 w 
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3?7 S \ dz dx dz 



(14) 

This is the equation for the dynamics of the film thickness - it is coupled to the particle 
volume fraction through the suspension viscosity i] s ((j)). 

Eq. (4) of can be recovered if we average ()14)) over the length of the entrainment region; 
neglecting the surface-tension terms for the moment and integrating (pPfj) from the pool 
(x = 0) to a distance L, we obtain 

dw pg_ 2 

— - axVu, = -a 2 — w (15) 

Ot Tj s 

where we have used 

w(x=0) — w(x=L) w(x=0) 



L L 

[cos9w 3 /3r] s } x=0 - [cos9w 3 /3r] s } x=L [cos9w 3 /3r], 



cm, 



,sja;=0 2 / 

-7 ~ ^ = Ot 2 W /T) s . 



As the film thickness rapidly decays to its asymptotic constant value (see Figs El and EJ) , it 
is sufficient to take L to be a few multiples of w(x=0). Eq. ()15|) is identical to Eq. (4) of 
our earlier paper j5j if the surface tension term in the latter is dropped. While the factors 
cki and o>2 are in general functions of z and t, they were treated as constants in first 
approximation: this is formally correct if we assume that the shape of the w(x) profile does 
not vary with z, even if its average does. As the most striking feature in the observations 
of was an axial variation in the particle concentration and the film thickness, this is 
a reasonable first approximation. Moreover, the results of this approximation are borne out 
by the full two-dimensional analysis in this work, as shown below. 

Surface tension was included incorrectly in [5|, making the naive guess that it entered simply 
at next order in gradients as awd 2 w/dz 2 . It is clear from (|14j) that it actually enters as 
aw 3 /3d A w/dz A , as acknowledged in our rebuttal to Nevertheless, the physical basis for 
the dynamical equation for the film thickness in is quite simple: it clearly follows from a 
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balance of gravity, shear forces and surface tension, and can hardly be described as "having 
no basis whatsoever" 

The equation governing is the mass balance for particles, 



d< t> 



(16) 



where j is the flux of particles due to shear-induced diffusion j^, Q|. The phenomena of 
shear-induced diffusion and migration have been addressed in several papers, and therefore 
does not require much description here. The shear-induced flux has two parts, one driven 
by a gradient in the concentration and the other by a gradient in the shear rate, 



j = -/ c a 2 7V0 - / s a 2 0V7- 



(17) 



where 7 is a measure of the shear rate, a is the particle size. The above can also be derived 

n 

by treating the particles and fluid as separate, but interacting, phases 15]; combining the 
balances of the particle phase mass and momentum yields an equation identical to (|17p to 
leading order in the gradients. As the flux depends only on the magnitude of the shear rate 
and not its direction, 7 must be a suitably chosen measure. If S is the rate of deformation 
tensor, the second invariant V2S : S is the natural choice for 7 jlj| in two dimensional flows. 
As we use the lubrication approximation, the only contributions to the shear rate rate come 
from S xy = S yx = dv x /dy. Further, since we assume no variation of the concentration in 
the y direction (across the thickness of the film), we use thickness-averaged shear-rate and 
velocity, 

rw pw 

v = l/w / vdy, (18) 
Jo 



7 = l/w / 'ydy, 
Jo 



in the mass balance. The particle mass balance then takes the form 



— + v x — + v z — = — 
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dt 



d_ 

dx 



, 2- d( t> , , 2 , °1 



dx 
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(19) 



Equations (|14|) and (|19|) together describe the dynamics of the suspension film. 
The parameters f c and f s are in general direction dependent - the diffusivity in the velocity 
gradient direction generally differs from that in the vorticity direction. However, we have 
taken them to be independent of direction in (|T7|) for simplicity. T hey are functions of the 
particle concentration. The function f c {4>) has been determined by 23| for migration in the 



vorticity direction in simple shear, but to our knowledge no direct measurement of f s exists. 
Using the suspension balance model [.15], the diffusivities may be related to the particle 
phase pressure, 

/ c = X2/9/(0)j/(0), f s = 2/9/(0)K0)/0 (20) 

where /(</>) is the hindered settling function in the sedimentation of particles in a suspension 
and p((f>) = 11/(777) is the dimensionless pressure function, 77 being the viscosity of pure fluid. 



The hindered settling function is taken to be f((f>) = (1 — </>) 5 [24(, and recent study of |2j 
gives p(4>) = rj s ((j))a((j)) with a{(f>) = 1.890 3 e 2 ' 34 ^. The suspension relative viscosity is taken 
as r] s = (1 -0/O.62)- 1 - 83 . 

We must emphasize that the above relations for f c and f s are only rough estimates, as there 
is some uncertainty in the data of j^j]; for instance, they have combined measurements of 
the stress in the suspension with an estimate of a normal stress in the particle phase to 
determine all the normal stresses. Due to the uncertaintly in their dependence on <ft, we 
treat f c and f s as constants in our model, equal to their respective values at = (fi p . The 
factor x i n the relation for f c is nominally equal to unity, but we have varied it to test the 
sensitivity of our predictions to the assumed forms above. 

A. Boundary conditions 

The solution in the entrainment region is matched with that of the static region when the 
film thickness in the former becomes very large B0Q. The matching is achieved by 
enforcing the boundary condition of that the curvature w" of the film must approach 
that given in (J3J) at w — > 00. In the solution of jlOj, the effect of gravity is neglected and the 
asymptotic behaviour for large w is indeed such that w" attains a constant value. However, 
an inspection of ()14j) shows that w'" — > 1/E as w — > 00, which gives a second derivative 
that varies linearly with x as w becomes large jll. 12 . 1^ |. An app roximate method has been 
used to satisfy the matching condition in some studies [3, but we use the following 
simple condition that has the same qualitative effect: we would like w'" to decay and w" 
to asymptotically approach the constant value given in (jSJ) as w — > 00, but since that isn't 
possible, we impose (J3J) at the point where w'" = 0. In other words, we patch the entrainment 
solution to the static solution where the curvature in the former is maximum. Consequently, 
the boundary conditions in dimensionless form at the lower boundary of the entrainment 
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region are 



w 



'2pg cos 9 



1/2 



w" = at x = 0. 



a 



(21) 



At the upper end of the entrainment region, the solution is matched with that in the constant 
film thickness region, and consequently the boundary conditions are 



w' = 0, w" = 0, vJ" = at x — ► oo. 



(22) 



For the particle concentration, we impose the condition that the concentration at the lower 
boundary is equal to that in the pool, i.e. 



at x = 



(23) 



and that the concentration gradient vanishes far from the pool, 



/ = at x — > oo 



(24) 



B. Base state solution 



The base state, whose stability we determine in the following section, is the steady state 
with no z- variation. Under these conditions, equations (jl4j) and (|19p reduce to 

dw d 



v 



w dx dx 
deb d 



p#cos# - a—s ) — 



d 3 w\ w 3 



dx dx 



f c a l— + f s a 



dx 



dx 



We now define the following dimensionless variables, 



w* = w/w s , x* = x/w s , z* = z/w s , r = v w t/w s 
v* = v/v w , 1 * = 1 w s /v w , r}* = ris(<fr)/Vs((j>p) 



and the dimensionless parameters 



e=—, £ 
wi 



a 



3v w r] s ((j) p 



1/2 



(25) 
(26) 



where w„ = \ j s sca l e for the film thickness. Henceforth, we work only in 

s \ pg cos v J ' J 

terms of these dimensionless variables, but drop the asterisks for the sake of convenience. 
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Equations (j23J)-(j2S|) then transform to 



(3u> V - „ W« t _ _ Em „X _ = _ (2?) 



£ (/cf ^ + /off + + ') 



7] 7] 

w 2 ' 

1 - (1 - Em'") — 



0, (28) 



7/ 

where the primes denote differentiation with respect to x, and M = d lnrj/dcf). The shear rate 
is \dv x /dy\ averaged over the film thickness, giving 7 = (1 — Hw'")3w/2. In dimensionless 
form, the boundary conditions (|2*T|) transforms to 



w 



" = (2/S) 1/2 , w'" = at x = 0. (29) 



and the others remain unaltered from ()22|) - (2~%). 
From (|2*7jl we recover the equation given by 
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12j for a homogeneous fluid by setting 0' 



to zero; further, if the gravitational body force is discarded, we recover the equation given 
by Q. 

The solution for the base state is obtained by integrating (|27j) - (j28j) numerically from x = to 
a suitably large distance x = Xoo. It is easily seen that the thickness and concentration vary 
as w(x) = w c + Ae~ kx , <f)(x) = <ft c + rAe~ kx as x — > 00, w c and C being the asymptotic film 
thickness and particle concentration, respectively. As a result, one may impose, at x = Xqo, 
the forms 

w = w c + a, w' = —ka, w" = k 2 a, w'" = —k 3 a, (30) 
= 0c + rex, 0' = —kra. 

The values of k and r are determined by substituting (p?0j) in (|27jl -(|28 jl and linearizing about 
a = 0. Upon eliminating r from the two equations, we get a quintic equation for k, and 
only one of its five roots is real and negative. The solution was obtained by integrating from 
x = Xoo to x = 0, with guesses for w c and C . The guesses were refined until boundary 
conditions (|2*T|) and ()23|) are satisfied. The solution is shown in Fig. El for a few values of S 
and e; it is clear that wq and 0o almost equal their asymptotic values a short distance from 
the suspension pool. 



12 



0.8 



0.6 



0.4 



0.2 







l\\ 


n 1 

U.l 




— i 


% 


10 















2 4 

x 



(a) 



0.4 



0.39 



0.38 





1 \ / 

\ \ / 

\ 

\ 

I 

- \ 

\ 
\ 

\ 1 


/ / 
/ / 
' / 
/ 

/ 

/ 

/ 
/ 

/ 

/ 


e 

— lO" 2 
lO" 3 

10 -< 


\ 1 
\ 1 
\ 1 

\ / 

\ / 







(b) 



FIG. 3: (a) The base state film thickness for three values of the dimensionless surface tension, (b) 
The base state particle concentration profile for three values of e with E = 1. The film thickness 
is insensitive to e in the range studied. In (a) e = 10 -2 and in (b) E = 1, and x = 1 i n both. 

III. LINEAR STABILITY 



We now allow perturbations to the base state so that the film thickness and particle con- 
centration have the form 



w = w (x) + w(x, Z,t), (f) = (j) (x) + z, t), 



(31) 



where Wq and <f>o are the base state solutions obtained in §TfB] Assuming the perturbations 
are small, i.e. w <C w , <C 4> , (fblj) and (fT^j) may be linearized about the base state. 
Finally, taking the plate to be infinitely long in the z direction, we may assume the solution 



to be in the form of normal modes in the z direction, w = w(x) e 
We then get for the linearized equations for the perturbations 
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with the shear rate perturbation 7 given by 
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k 2 D )w - fl - Ew, 
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(34) 



where .D represents the operator d/dx, and c = u> c — w^. 

Boundary conditions for f!32l) and (f3*3*|) are obtained by linearizing ()21j1 - ()24j) , 



Dw = D 2 w = D 3 w 



at x = 
= at a; = 



X r 



(35) 
(36) 



Before solving the eigenvalue problem (JH2J)-((HEJ), it is instructive to consider the simplest case 
of a film of uniform thickness dragged by (or draining from) a flat plate, as is the situation far 
above the pool. As this problem is translation-invariant in x, we may Fourier-transform in 



the x direction as well as z, w 



w e 



ik x x+ik z z+st 



, e ik x x+ik z z+st _ fj^g ii near i zec i equations 



describing this simplified problem are 



sw 



{-Zwl(k 2 x + k 2 z ) 2 + i (-1 + 3wl) k x ) w - {iMw'lK} } 
hw + iY>w 2 k x {k 2 x + k\) - tik 2 x + k 2 z )^^ [l + iEw k x (k 2 x + k 2 ) 



-'/. (p 4- b 2 \ < ' 2 ' 2\3fsVU 

i{\ - wl)k x + e{k 2 x + k* 



w 



(37) 
(38) 



3w 



The eigenvalues of the above equations were determined and expanded as a series in k x and 
k z using the MATHEMATICA software package. The results up to second order in the wave 
numbers are 



T A/" 1 / 2 ?/; 2 h 1 / 2 + W ° 
T-/V W Q K X =■= 4 jyl/2 



kl /2 -- A ew {f c -UMs) (k 2 x + k 2 z ) 



+i 



±M l ' 2 wlkl l2 + 



•1 + 2w 2 ) k x ± 



(39) 



4AA 1/2 

When k x = 0, the result is essentially what was given in there is an instability when 
the discriminant E = f c — A/"0o/s turns negative, with the instability growth rate varying as 



—Ek 2 . The mechanism of this instability, as given in [5], is the following: a rise in <ft at some 
axial position causes an increase in the local viscosity rj s (4>), leading to a fall in the local 
shear-rate there, in part by increasing the local film thickness w. Shear-induced diffusion 
causes particles to migrate towards regions of lower 7, further increasing in the region and 
therefore driving an instability. We note here that [16j analyzed the linear stability of wall- 
bounded plane Couette flow, and predicted that the system is stable if driven at constant 



14 



shear rate, but unstable at constant shear stress provided —f c /7 < E < 0; the condition of 
instability for a uniform film is clearly less stringent. 

There is also an instability at non-zero k x , which varies as k]/ 2 for small k x , and is present 
even when k z = 0. This "pile-up" instability, is not driven by shear-induced migration and 
hence does not involve the discriminant E; it arises from the concentration dependence of the 
viscosity and the advection of concentration perturbations in the x direction. The stabilizing 
effect of surface tension comes in at fourth order in the wavenumbers in all directions. 
When the region near the pool is considered, the problem is not translation-invariant in 
x, and we therefore return to the full problem l[32j) - l|H6|) . The eigenvalue problem is solved 
numerically by a Chebyshev collocation spectral method. The physical domain x = (0, a^) 
is mapped to the computational domain £=(—1,1) by the transformation 

x = I -. (40) 

The linearized equations are enforced at the N collocation points, which are the roots of 
the N th order Chebyshev polynomial T/v(£) = cos(iV cos -1 £). The parameter I is used to 
increase the density of collocation points near the pool, where there is rapid variation of the 
variables. The variables w and <p are expanded in Chebyshev polynomials, w = Y$Li a iTi(£), 
<p = Y,iLi biTi(£), reducing to a set of 2N linear algebraic equations for the a, and 

6j. The eigenvalues of the coefficient matrix give the growth rates s r of the perturbations. 
The real part of the leading eigenvalue s l r is termed the leading growth rate, and its variation 
with k z and the dimensionless parameters in the problem are presented in Figs 0]- E] below. 
For all our computations, was at least 16, I = 10 and N = 81, which yielded eigenvalues 
with accuracy greater than 10~ 5 . The results presented below are for <p p = 0.4; our results 
for lower <ft p are qualitatively the same, but the numerical accuracy of the eigenvalues is 
reduced. 

Figs E] and El show the dependence of s l r on the parameters £ and x- Fig. |U1 shows the profile 
of the eigenf unctions for <ft and w as functions of x for k z — 1, showing that they are localised 
close to the pool. Several important points emerge from the figures: (i) Instability is always 
present, irrespective of the value of the dimensionless surface tension E. This is consistent 
with the discussion and results for the straight film . (ii) The instability persists down to 
k z = because, as seen in Fig. El variations along x are always present in the eigenfunctions 
so that the pile-up instability always arises; the mechanism for instabilities with appreciable 
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FIG. 4: The leading growth rate s r as a function of the wavenumber k z for some representative 
values of the dimensionless surface tension E. Here the particle concentration in the pool is 4> p = 0.4 
and the parameter x is 0.5, yielding a value of 4.47 for discriminant E, and e = 10~ 2 . 

variation along z is, as we have said above, essentially that outlined in riding on the pile- 
up instability. If the discriminant E is below a critical value E c , the growth rate above the 
baseline value at k z = varies as k\ for small k z , as stated in |5]. While the one-dimensional 
analysis of \^ predicted E c = 0, our full two-dimensional here predicts E c > (see Fig. EJ). 
(hi) Surface tension tends to stabilize modes at large wavenumber and hence determines the 
range of unstable wavenumbers. This can be seen clearly in Fig. JU (iv) Changing the value 
of the discriminant E = (f c — N<fiofs) has the effect discussed in |5[: a large positive value of 
that parameter, controlled here by the parameter x, greatly reduces the growth rate (which 
would in fact go negative if it were not for the effect of variation along x). This can be seen in 
Fig. El where the curve for \ = 0.45 corresponds most closely to a case dominated by modes 
varying along z, and hence shows the increase in growth-rate with increasing wavenumber. 
(v) For typical parameter values considered, the most unstable wavenumber is at k z ~ 1, 
which implies a wavelength of ~ 2ttw s . For the experiments of 0, w s in the range 1-3 mm, 
and hence the wavelength of the fastest growing perturbations would be around 1 cm, which 
is within their observed range of the distance between adjacent bands. 
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FIG. 5: The leading growth rate s r as a function of k z for different values of the factor x ( see 
Eg. I20|). The values of the discriminant E = (f c — M4>of s ) for the four cases are 57.76, 31.11, 4.47 
and —0.86 (in the order of decreasing \)- Parameter values are (j) p = 0.4, S = 1, and e = 10 -2 . 

IV. SUMMARY AND CONCLUSION 

We have shown that a film of suspension coating a wall rising from a suspension pool can 
display an instability to perturbations that cause segregation of particles in axial bands. 
This instability is particularly pronounced when the discriminant E, which depends on the 
dimensionless shear-induced diffusivities and the suspension viscosity, is below a critical 
value E c . This bears out the results we obtained in an earlier paper j5j, now in a complete 
analysis in which the equation governing the dynamics of the film thickness is derived from 
the Stokes equations and the kinematic condition at the free interface, and a fully two- 
dimensional stability analysis carried out. The instability comes from the coupled dynamics 
of film thickness w and particle concentration cf). The condition for instability we predict is 
less stringent than the condition —f c /7 < E < for plane Couette flow at constant shear 
stress given by jig . 




We observe two kinds of instabilities: the first is the axial instability proposed in our ear- 
lier one-dimensional treatment, in which axial variations of the film thickness and particle 
concentration are amplified with time. The instability growth rate for this instability varies 
as the square of the axial wavenumber k z for small k z , precisely as in the one-dimensional 
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FIG. 6: The eigenfunctions for the particle concentration, <f>(x), and the film thickness, w(x), for 
k z = 1. Notice that they are restricted to a fairly small range of x. Parameter values are 4> p = 0.4, 
S = 1, e = 1(T 2 and % = 1- 

model of 0. The mechanism of this instability is that a rise in <fi at some axial position 
causes an increase in the local viscosity r] s ((p), leading to a fall in the local shear-rate 7, in 
part by increasing w; shear-induced diffusion causes particles to migrate towards regions of 
lower 7, further increasing <fi in the region and therefore driving an instability. In the second 
"pile-up" instability, variations along the x direction (i.e. of wall motion) are amplified, and 
the growth rate is finite even at k z = 0. This mechanism of this instability is that a rise in 
cf) at one position results in a rise in the local viscosity and hence a reduction in the veolcity. 
This reduces the advection flux of particles from that spot, causing a further accumulation 
of particles downstream. 

While we find that the eigenfunctions always have a significant x-dependence, and therefore 
the instability is inherently two-dimensional, the dependence of the instability growth rate 
on the axial wavenumber is as in the one-dimensional model j^jwhen E < E c . Surface tension 
plays a stabilizing role at larger wavenumbers, but in a manner different in detail from that 
proposed in Q] . To see whether this mechanism is indeed responsible for the observations 
of |ll [J, independent measurements of the dimensionless shear-induced diffusivities f c and 
f s (cf. Eq. |T7j) are required. 
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LIST OF FIGURE CAPTIONS 



Figure 1: Schematic diagram of film coating in a horizontal rotating cylinder. 

Figure 2: Coating of a vertical plate rising with velocity v w from a pool of suspension. 

Figure 3: The base state film thickness for three values of the dimensionless surface tension, 
(b) The base state particle concentration profile for three values of e with X = 1. The film 
thickness is insensitive to e in the range studied. In (a) e = 1CT 2 and in (b) £ = 1, and 
X = 1 in both. 

Figure 4: The leading growth rate s r as a function of the wavenumber k z for some represen- 
tative values of the dimensionless surface tension S. Here the particle concentration in the 
pool is 4> p = 0.4 and the parameter x is 0.5, yielding a value of 4.47 for discriminant E, and 
e = 10- 2 . 

Figure 5: The leading growth rate s l r as a function of k z for different values of the factor x 
(see Eq. I2UJ) . The values of the discriminant E = (f c — Af(f>of s ) for the four cases are 57.76, 
31.11, 4.47 and —0.86 (in the order of decreasing x). Parameter values are <p p = 0.4, E = 1, 
and e = 10~ 2 . 

Figure 6: The eigenfunctions for the particle concentration, 4>(x), and the film thickness, 
w(x), for k z = 1. Notice that they are restricted to a fairly small range of x. Parameter 
values are (fi p = 0.4, S = 1, e = 10~ 2 and % — 1- 
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